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Abstract — The differentiation by integration method with Ja- 
cobi polynomials was originally introduced by Mboup, Join and 
Fliess 1 22 1, |23|. This paper generalizes this method from the 
integer order to the fractional order for estimating the fractional 
order derivatives of noisy signals. The proposed fractional order 
differentiator is deduced from the Jacobi orthogonal polynomial 
filter and the Riemann-Liouville fractional order derivative 
definition. Exact and simple formula for this differentiator is 
given where an integral formula involving Jacobi polynomials and 
the noisy signal is used without complex mathematical deduction. 
Hence, it can be used both for continuous-time and discrete- 
time models. The comparison between our differentiator and 
the recently introduced digital fractional order Savitzky-Golay 
differentiator is given in numerical simulations so as to show its 
accuracy and robustness with respect to corrupting noises. 

I. INTRODUCTION 

Fractional models arise in many practical situations (Q, 
11331 for example). Such fractional order systems may also be 
used for control purposes: CRONE control is known to have 
good robustness properties (see lf26ll . 1271 . J28], ED). In order 
to implement such controller one needs to have a good digital 
fractional order differentiator from noisy signals, which is the 
scope of this paper. 

The fractional derivative has a long history and are now 
very useful in science, engineering and finance (see, e.g., 1251 . 
ifTTl . IE)). The fractional order differentiator is concerned 
with estimating the fractional order derivatives of an unknown 
signal from its noisy observed data. Because of its importance, 
various methods have been developed during the last years. 
They are divided into two kinds: continuous-time model (see, 
e.g., |32), 12, (29l) and discrete-time model (see, e.g., |38l , 
1351 , ||39l ). Nevertheless, the real case of signals with noises 
was somewhat overlooked. In order to resolve this problem, 
an optimization formulation using genetic algorithms was 
proposed in 112011 to reduce the noise effect in the estimations 
of the fractional order derivatives. But, the complex math- 
ematical deduction restricts its application. A novel Digital 
Fractional Order Savitzky-Golay Differentiator (DFOSGD) 
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was introduced in [3|, which was deduced from the Riemann- 
Liouville fractional order derivative definition iPTO (p. 63) and 
the Savitzky-Golay filter |36L ED. Using the Savitzky-Golay 
filter guarantees its accuracy and simplicity for estimating 
the fractional order derivatives of noisy signals. However, let 
us recall that the Legendre polynomial filter J30i] is more 
recommended than the Savitzky-Golay filter for reasons of 
simplicity and speed. In particular, it has advantages of being 
suitable for irregularly spaced or missing data. 

The method of differentiation by integration uses an integral 
of an unknown noisy signal to estimate the integer order 
derivatives of this signal. For free of noise signals, this method 
was firstly studied by Cioranescu J4) (1938) and became 
well known for the Lanczos generalized derivative fPJl (p. 
324) (1956). The Lanczos generalized derivative proposed 
an integral of a noisy signal and the Legendre orthogonal 
polynomial to estimate the first order derivative of the signal. 
Recently, in the noisy case, Mboup, Join and Fliess 11221 . lf23l 
applied an algebraic setting to estimate high order derivatives 
by integration, where Jacobi orthogonal polynomials were 
introduced in the integral. Hence, we call the obtained differ- 
entiator Jacobi differentiator. Moreover, it was shown in lfl5l 
that the Jacobi differentiator greatly improved the convergence 
rate of the Lanczos generalized derivative. Very recently, it 
was shown in IfTTl that the Jacobi differentiator could also 
be obtained by taking the derivative of the Jacobi orthogonal 
polynomial filter considered as the extension of the Legendre 
polynomial filter 1301. l24l. 

Let us recall that the algebraic differentiation method used 
in 1221 . 1231 was introduced in J9] and also analyzed in |[T4l . 
lfT31 . 02), IfTTl . 1341 . where the used algebraic manipulations 
were inspired by the algebraic parametric estimation methods 
ifTOl . OTi . lfT8ll . Ifl9l . Additional theoretical foundations can 
be found in (6], |8). In particular, by using non standard 
analysis techniques, Fliess J6), showed that these methods 
exhibit good robustness properties with respect to corrupting 
noises without the need of knowing their statistical properties. 
However, these methods have not been used to estimate 
fractional order derivatives. 

The aim of this paper is to generalize the differentiation 
by integration method from the integer order to the fractional 
order for estimating fractional order derivatives. In Section 
UTI we recall the method of integer order differentiation 
by integration with Jacobi polynomials. Then, we deduce 
a fractional order differentiator from the Riemann-Liouville 
fractional order derivative definition and the Jacobi orthogonal 
polynomial filter. This differentiator is exactly given by an 



integral of Jacobi polynomials. Hence, we call it fractional 
Jacobi differentiator. In Section [TTTJ we compare the fractional 
Jacobi differentiator to the DFOSGD in some numerical sim- 
ulations. Finally, we give some conclusions and perspectives 
for our future work in Section |IV] 

II. METHODOLOGY 

Let y s = y + vo be a noisy signal observed in an interval 
I = la, b] C R of length h = b - a, where y G C(I) and the 
noisqj w is bounded and integrable. We are going to estimate 
the a th (a G R+) order derivative of y by using its observation 

A. Differentiation by integration with Jacobi polynomials 

In this subsection, we recall the method of integer order 
differentiation by integration with Jacobi polynomials. This 
method was introduced in ll22l . 1231 and studied in lfT4l . lfT31l . 

ED, EO, M- 

The n th (n G M) order Jacobi orthogonal polynomial 
defined on [0, 1] is given as follows (see E p. 775) 



•s—^ I n + /i \ / n + K 
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where /i, k g] — l,+oo[. Let / and g be two functions 
belonging to C([0, 1]), then we define the scalar product 
('' ')u k °f m ese functions by (see HI p. 774) 



</(•), ff(0>u,„= / ^, K (t)f(t)g(t)dt, 



(2) 



where u> MiK (t) = (1 — t)^t K is the associated weighted 
function. Thus, by denoting its associated norm by || • \\^, K , 
we obtain 
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+ k + n + 1) T(n + 1) (2n + /i + k + 1) ' 

where T(-) is the classical Gamma function (see E p. 255). 

Let us ignore the noise for a moment. Then, we define an 
N th (N G N) order approximation polynomial of y by taking 
its truncated Jacobi orthogonal series expansion 
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If we take k = /i = 0, then the Jacobi orthogonal polynomials 
become the Legendre orthogonal polynomials. Hence, this Ja- 
cobi polynomial filter l24l can be considered as the extension 
of the Legendre polynomial filter l30l . 

'More generally, the noise is a stochastic process, which is bounded with 
certain probability and integrable in the sense of convergence in mean square 
(see Q21). 



If y G C n {I) with n G N, then we take the n th order 
derivative of the polynomial D£' K N y(a + h-) as an estimate 
of the n th order derivative of y lfl7l : Vi G [0, 1], 

{DZ, K , N y(a + ht)} (5) 
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For any t G [0,1], this differentiator can be expressed as 
follows Q7) 

D%l, KtN y(a + ht) = 
1 f 1 

Q,j,, K ,n,N(T, t) y(a + hr) dr, 
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Finally, we replace y in (O by y" 5 . Consequently, the 



order derivative of y can be estimated by an integral of Jacobi 
polynomials. The noise effect on obtained estimations was 
analyzed in E), ED, ED, E3, EQ. In the next subsection, 
we are going to generalize this differentiation by integration 
method from the integer order to the fractional order. 

B. Fractional order differentiation by integration 

Similarly to 0, we are going to take the Riemann-Liouville 
fractional order derivative of our approximation polynomial 
so as to get a fractional order differentiator. The Riemann- 
Liouville fractional order derivative (see BP p. 62) is defined 
as follows: M x G Rl, 



1 



(x-t) 1 -*- 1 f(r)dT, (7) 



T(l — a) dx l J 

where < I - 1 < a < I with / G N*. 

From now on, we denote the a fh (a G R+) order derivative 
of /by 



d° 



V* G R;, D«f(x) = — {f(x)} = fM(x) 



(8) 



Hence, if we take f(x) = x n with n G N and x G R5_, then 
we obtain (see 1311 p. 72) 



^ rn,_ r(n+l) 

dx a x ' T(n + 1 - a) 
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By using (O, we obtain the following lemma. 

Lemma 1: The a th (a G R+) order derivative of the n th 
order Jacobi orthogonal polynomial P^ 1 '^ defined in El is 
given as follows: Vt G]0, 1], 

£ {*'"'(')} = 

" V(n 1 A- \\ 

/ j / j C IJ;K,n,j,l - - - n C 
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where = (-l)^f) (V)- 

Proof. By applying the binomial theorem to (Q3, we get 

n n—j 

where c^,^ = (-l)'^") (V)- Hence > this P roof 
can be completed by using (|9]l and the linearity of the 
fractional order differentiation (see ||3T1 p. 91). □ 

Then, we give the following proposition. 

Proposition 1: Let y s = y + w be a noisy observation of 
y in an interval / = [a,b] C K of length h = b — a, where 
y 6 C(I) and the noise w is bounded and integrable. If the 
a th (a £ R+) order derivative of y exists, then a fractional 
order differentiator, called fractional Jacobi differentiator, is 
defined as follows: Vt e]0, 1] 



1 f 1 

J Qn,K, a ,N{.T,i)y 5 {a + hT)dT, 
where ft e i;, JV e N, ft k e] - 1, +oo[, 
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with c^,, = (-lyCTOG^HV) and HP^H^ is 
given in ©. 

Proof. Let us take the a th order derivative of the polynomial 
Dfr' K N y(a + h-) defined in @. By using the scale change 



property of the fractional order differentiation (see 1125] p. 76) 
we obtain: Vt e]0, 1], 

The linearity of the fractional order differentiation gives us 
that 



.a / t 
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Finally, this proof can be completed by using Lemma Q] and 
substituting y by y 5 in (fT5l l. □ 
Consequently, by calculating the integral of the noisy signal 
y 5 and a sum of the Jacobi polynomials we can estimate the 
value of t/ q ' at each point a + ht in the interval ]a, b] for 



each t g]0, 1]. This integral corresponds to a convolution in 
the discrete case. 

If we fix the value of t to 1, then the fractional Jacobi 
differentiator D^f 1 K N y s (a + ht) only estimates the value 



of yW at point b. Thus, if we increase the length of the 
interval I = [a,b] C K, then we can estimate the other values 
of y^ a '. Consequently, the fractional Jacobi differentiator can 
also be considered as a pointwise causal differentiator which 
is useful in many on-line applications. Moreover, we can fix 
the parameters k, fj, and N such that the sum of the Jacobi 
polynomials can be explicitly calculated by off-line work. 
Indeed, the fractional Jacobi differentiator only needs a simple 
convolution for on-line applications. The computation time is 
significatively improved. 

III. SIMULATIONS RESULTS 
The comparisons between the DFOSGD and some existing 
fractional order differentiators were shown in [3], where the 
DFOSGD was better than the others. In this section, by taking 
one numerical example considered in O we compare the 
fractional Jacobi differentiator to the DFOSGD so as to show 
its accuracy and robustness with respect to corrupting noises. 

A. Numerical case 

From now on, we assume that y s (xi) = y(xi) + cw(xi) is 
a sequence of uniformly sampled noisy data of y with Xi = 
a + T s i and T s = jj for i = 0, ■ • • , M. The noise cw(xi) is 
assumed to a zero-mean white Gaussian iid sequence, and c 6 
R+. In this discrete case, we apply the trapezoidal numerical 
integration method in the fractional Jacobi differentiator to 
approximate the integral. Then, let us recall that the values of 
y {a \xi) for % = 0, • ■ ■ , M are estimated by the DFOSGD as 
follows 

1 , 
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Hence, the idea is to obtain an N th order approximation 



/ 



polynomial by using Savitzky-Golay algorithm j36l , ll37l 
from a subsequence Yg, and then to calculate the a th order 
derivative of this polynomial at each point Xi. 
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Fig. 1. Signal y and noisy signal y . 



Fig. 2. Absolut estimation errors in noise-free case. 



B. Example 

In this example, we assume that y s (xi) — sm(5xi)+cw(xi) 
where Xi e I = [0, 4] for i = 0, • ■ • , 1000 = M, and c = 0.25 
is adjusted in such a way that the signal-to-noise ratio SNR = 

101 °gio ( ^\cL%% ) is ec i ual to s-W-ft = 10dB - The nois y 

signal is shown in Figure Q] 

Firstly, we estimate the half order derivative of y. For 
the DFOSGD, we take N = 14 and 9 = 5 which are 
the same values as the ones used in Q. For the fractional 
Jacobi differentiator, we take N = 14 and k = /i = 0. The 
associated absolute estimations errors in the noise-free case 
and in the noisy case are given in Figure [2] and Figure [3] 
respectively, where the black dash-dotted lines represent the 
errors obtained by the fractional Jacobi differentiator and the 
red dotted lines represent the ones obtained by the DFOSGD. 
Moreover, we can see in Figure|4]the associated absolute noise 
error contributions. Secondly, we estimate the derivative j/ Q ) 
with a = 1.5 by using the DFOSGD and the fractional Jacobi 
differentiator with N = 14, 9 = 5 and k = /x = 0. The 
obtained estimation errors are shown in Figure |5j Figure [6] 
and Figure [7] Consequently, we can observe that the fractional 
Jacobi differentiator is more accurate and more robust with 
noises than the DFOSGD in the noise-free case and in the 
noisy case. 

IV. CONCLUSION 

In this paper, we propose a fractional Jacobi differentiator 
which is a fractional order differentiator exactly given by 
an integral formula. This differentiator is deduced from the 
Jacobi orthogonal polynomial filter and the Riemann-Liouville 
fractional order derivative definition. It can accurately and 
easily estimate the fractional order derivatives of noisy signals. 
There are some parameters h, fx, k and TV on which the 
fractional Jacobi differentiator depends. We do not give any 
analysis on the influence of these parameters on the estimation 
errors. However, similarly to the integer order differentiation 



Absolute Estimation errors 



0.5 
0.45 

0.4 
0.35 

0.3 
0.25 

0.2 
0.15 

0.1 
0.05 




I 
I 

' I 
' I 
I I 

i ' 

i - 

i 



■ DFOSGD 1 

D (c.0.5)„S I 



\ ■ 







0.5 



1.5 



2 
Time 



2.5 



3.5 



Fig. 3. Absolut estimation errors in noisy case. 
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Fig. 4. Absolut noise error contributions. 
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Fig. 5. Absolut estimation errors in noise-free case. 
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Fig. 6. Absolut estimation errors in noisy case. 
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Fig. 7. Absolut noise error contributions. 



by integration case (see lfl6l . [fT71 . fl31 . lfl4l ). this study can 
be easily carried out in a further work. Moreover, we do not 
use the causal propriety of the fractional Jacobi differentiator 
in the numerical simulations. It will be used in some interest- 
ing applications where the sampling data may be irregularly 
spaced. 
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